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ABSTRACT 

We investigate the observational signatures associated with one of the proposed for- 
mation scenario for the recently discovered highly eccentric binary millisecond pulsar 
(MSP) PSR J1903 + 0327 in the galactic plane. The scenario requires that the MSP 
to be part of a hierarchical triple (HT), consisting of inner and outer binaries, expe- 
riencing the Kozai resonance. Numerical modeling of a bound point mass HT, while 
incorporating the effects due to the quadrupolar interactions between the binary or- 
bits and dominant contributions to the general relativistic periastron precession in 
the inner binary, reveals that, at the present epoch, the orbital eccentricity of the bi- 
nary MSP should decrease for reasonable ranges in the HT parameters. The estimated 
decrements in the orbital eccentricity of the inner binary are few parts in 10^, substan- 
tially higher than the reported accuracies in the estimation of the orbital eccentricity 
of the binary MSP, while employing various general relativistic timing models for iso- 
lated binary pulsars. For wide ranges in the allowed orbital parameters, the estimated 
rate of change in the eccentricity of the inner binary is orders of magnitude higher 
than the value recently measured by the pulsar timing analysis. Therefore, we rule 
out the scenario that the MSP is part of a HT undergoing the Kozai oscillations. The 
origin of this system in a typical globular cluster is also shown to be less likely than 
inferred in the discovery paper. 

Key words: stars: neutron — pulsars: individual (PSR J1903+0327) - methods: 
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1 INTRODUCTION 

The recent discovery of a MSP PSR J1903-I-0327 in a highly 
eccentric orbit around a solar mass companion in the galac- 
tic pla ne forces one to explore unconventional formation sce- 
narios (|Champion et al. 20081 ) . This is because it is very un- 
usual for a galactic binary pulsar to have spin periods shorter 
than 10 ms and a large orbital eccentricity. In the standard 
model for the formation of the galactic MSPs, accretion of 
mass and angular momentum from a binary companion is 
required for the spi n-up of the pulsar , necessary to create 
short spin periods (|Alpar et al. 1982 ). Due to the strong 
tidal effects operating in the binary during the accretion 
phase, MSP bi naries are exp ected to be in almost perfect 
circular orbits l|Phinnev 19921 ). 

Accurate radio timing of the MSP reveals that the ec- 
centricity of the binary, having an orbital period of around 
95 days, to be ^ 0.44 and the reported accuracy of its esti- 



mation is ~ 5 parts in lO" , while using a theory- independent 
general relativistic binary pulsar timing model. Further, the 
estimated mass of the MSP is 1.74 ± 0.04Mq. Infrared ob- 
servations with Gemini North telescope yielded a possible 
main-sequence companion star for the MSP. The presence 
of the main-sequence star in the MSP field prompted its 
discoverers to suggest triple star interactions to explain the 
unusually high orbital eccentricity of the binary MSP. Let us 
briefly summarize other plausible formation scenarios for the 
binary system. Due to its distinct spin and spin-down rates. 
Champion et al. (2008) have almost discarded the "born- 
fast" scenario for this pulsar and strongly favored a binary 
origin for it. They have also speculated a globular cluster 
(GC) origin for the system, involving the ejection of the bi- 
nary MSP from a GC or the disruption of the cluster itself. 
These cases are also invoked to provide some explanation 
for the high orbital eccentricity of the system. 

One of the most plausible formation scenario for 
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PSR J1903 + 0327 requires a HT configuration involving the 
MSP, its binary companion, likely to be a massive white- 
dwarf, and the main-seq uence star, observed in the infrared 
in the field of the MSP (|Champion et al. 200 j ). In this sce- 
nario, the MSP is part of a primordial HT system com- 
posed of the pulsar and the massive white-dwarf. The main- 
sequence companion star is in a much wider and highly in- 
clined orbit around the inner binary. The tidal interactions 
between the inner and outer binaries produce eccentricity 
oscilla tions in the inner binary through the Kozai mecha- 
nism (|Kozai 19621 ). Champion et al. (2008) estimate that 
it is reasonable to expect around two HT systems among 
the known MSP binaries in the galactic disk, based on an 
estimate that gave roughly 4% chance for a given neutron 
star-white dwarf binary to be part of a triple. With the help 
of Ford, Kozinsky & Rasio (2000), the authors also noted 
that it is possible for the inner binary of a highly inclined 
HT to have eccentricity ~ 0.44 during 20% of its Kozai os- 
cillation period, when its masses are 1.74, 1.05 and 0.9 M©, 
respectively. These estimates make a HT configuration the 
most plausible formation scenario for PSR J1903 -I- 0327. 

In this article, we explore the physical consequences of 
Kozai oscillations operating in the system and its relevance 
for the on-going timing of PSR J 1903 -1-0327. Using reported 
values for the orbital eccentricity and the longitude of pe- 
riastron, extractable from the Table. I in Champion et al. 
(2008), and reasonable estimates for various parameters nec- 
essary to define a HT, we demonstrate that, at the present 
epoch, eccentricity of the binary MSP should decrease. This 
is achieved by modeling the system to be part of a bound 
point mass HT, while incorporating the effects due to the 
quadrupolar interactions between the binary orbits and the 
dominant order general relativistic periastron precession in 
the inner binary. Typically, the accumulated eccentricity de- 
crease in an year for the inner binary is substantially higher 
than the reported accuracy with which the timing eccentric- 
ity is estimated (few parts in 10^ against few parts in 10^). 
Further, we note that the recently measured rate of change of 
the orbital eccentricity ~ 10~^^s~^ (Freire, priv. commun.) 
unequivocally rule out the above described HT scenario for 
the MSP binary. We also make an estimate for the spin pe- 
riod derivative, assuming that the MSP is in a triple system 
and find it to be inconsistent with observations. 

Our investigations pertaining to other formation scenar- 
ios outlined in Champion et al. (2008) reveal that the direct 
observation of the massive white-dwarf companion through 
its optical band radiation may not be possible in the near fu- 
ture. Additionally, the origin of this binary MSP in GCs may 
be subject to more restrictions than has been realized so far. 
This is because the MSP system like PSR J 1903 -I- 0327 can 
originate only from GCs having special properties, thereby 
reducing the overall probability of this pulsar's origin in the 
galactic GC systems. 



2 ESSENTIAL MODEL TO DESCRIBE THE 
DYNAMICS OF HIERARCHICAL TRIPLES 
EXPERIENCING THE KOZAI 
OSCILLATIONS 

The HT formation scenario for PSR J 1903 + 0327, suggested 
by its discoverers, allows us to treat the system to consist 



of two binaries in quasi-Keplerian orbits that are highly in- 
clined to each other. The inner binary contains the MSP and 
its stellar-mass companion, possibly a massive white-dwarf, 
with masses mo and mi. The stellar-mass main-sequence 
star having mass m2 forms an outer binary with the center 
of mass of the inner binary. Let us denote the eccentricities, 
semi-major axes and the arguments of the periastron (with 
respect to their lines of nodes) of the inner and outer bi- 
naries by 61,62,01,02,51 and 32, respectively and let i be 
the mutual inclination angle between the two orbits. We 
probe the temporal evolution of the inner binary with the 
help of secular perturbation theory, applicable to Newtonian 
HTs containing point masses, while including the dominant 
quadrupolar order interactions between the two orbits. In 
other words, the dynamical equations that we invoke are 
accurate to order (01/02)'^, where (01/02) is the small pa- 
rameter in the perturbative expansion. We also incorporate, 
in an ad-hoc manner, the dominant order general relativis- 
tic efi'ect that causes the periastron of an isolated compact 
binary to advance through the dynamical equation for the 
argument of periastron of the inner binary. The relevant 
equations providing secular temporal evolution for the ec- 
centricity and the argument of periastron of the inner binary 
(|Blaes et al. 20021 ') read 
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where Mi — (mo -I- mi), the total mass of the inner binary 
and 9 = cos i. The quantity C2 and the magnitudes of the 
angular momenta Gi and G2 of the inner and outer binaries 
are given by 
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where rji is the symmetric mass ratio of the inner binary, 
given by rji = momi/Mf. Note that the Newtonian con- 
tributions to Eqs. ([1} originate from certain 'doubly av- 
eraged' Hamiltonian, which is derivable from the usual 
Ha miltonian for a H T at the quadrupolar interaction or- 
der (|Ford et al. 20001 ). The 'doubly averaged' Hamiltonian, 
suitable for describing secular (long-term) temporal evolu- 
tion of a HT, is independent of the mean anomalies of the 
inner and outer orbits. This implies that their respective 
conjugate momenta and hence the semi-major axes, oi and 
02, are constants of motion. Moreover, the justifiable neglect 
of the radiative losses to the orbital energy and angular mo- 
mentum of the inner binary due to the emission of gravita- 
tional waves implies that there are no reactive contributions 
to dai/dt and dei/dt. 

When one neglects general relativistic contributions and 
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let G2 » Gi in Eqs. ffl, the resulting equations allow an 
analytic solution (jKozai 19621 ). It is also possible to con- 
struct an approximate integral of motion in terms of ei and 
6 that allows one to classify the dynamical behavior of a HT 
with the help of trajectories in the the phase space defined 
by ei and cos gi. It turns out that if the mutual inclination 
angle is fairly high and in a certain window, the orbital ec- 
centricities experience periodic oscillations over time-scales 
that are extremely large compared to the respective orbital 
periods. The above effect, usually referred to as the Kozai 
resonance, arises due to the tidal torquing between the two 
orbits. The Kozai resonance can force initially tiny eccen- 
tricity of the inner binary to oscillate through a maximum 
value, given by ei™^'' ~ (l — | cos^io)^''^ where io is the 
initial value for the mutual inclination angle. Due to the ob- 
vious restriction, namely | cos ^o < (3/5)"^^'^ |, io is required 
to lie in the range 39° - 141° (|Blaes et al. 2002h . 

The general relativistic periastron advance of the inner 
binary, in principle, can interfere with the Kozai resonance 
and even terminate the eccentricity oscillations. This is be- 
cause the extra contributions to dgi /dt can indirectly affect 
the evolution of ei. The following useful criterion, derived 
in Blaes, Lee & Socrates 2002, may be used to infer the 
possibility for this not to happen. 
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where Mi is Mi/MQ,ai = oi/Lq and Lq = 1.476625 km. 
The above inequality is obtained by equating the right hand 
side of Eq. (|la|) to zero, after neglecting the much smaller 
C2/G2 contributions to dgi/dt and demanding that the re- 
sulting expression for cos^ i remains positive. 

Another constraint for a2/a\ can be obtained 
by invoking an emph erical relation obtained by 
l|Mardling fc Aarseth 20011 ). relevant while discussing 
the stability of Newtonian coplanar prograde or- 
bits in HT configuration s. The empherical criterion of 
iMardling fc Aarseth 200ll ) reads 
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Following Blaes, Lee & Socrates (2002), we treat the above 
inequality to be rather conservative as the inclined orbits, 
relevant for our investigation, are expected to more stable 
than the coplanar triples of Eq. Q . 

An initial estimate for 62, appearing in the above in- 
equalities, can be obtained by equating the general rela- 
tivistic periastron precession timescale of the inner binary 
to the characteristic time scale for the Kozai oscillations. 
The period for general relativistic periastron precession is 
deductible from the reported period for advance of perias- 
tron ~ 10" years, which assumes that the MSP binary is 
isolated. An approximate expression for the period of the 
Kozai oscillation reads 
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where Pj is the orbi tal period of the inner binary 
l|Mazeh fc Shaham 19791 ). With these inputs, for a HT sys- 
tem with Mi ~ 2.79MQ,ai ~ 0.211 A.U, 02 ~ 600 x 
ai,m2 ~ O.9M0 and Pi ~ 0.261yr, the likely values for a 



HT listed in Champion et al. (2008), the above prescrip- 
tion gives 62 ~ 0.979. Not surprisingly, the above value for 
62 fails to satisfy both the inequalities. However, we observe 
that a slightly higher value for 62, say 62 ~ 0.984, makes sure 
that the above two inequalities are satisfied. This is how we 
prescribe a suitable value for 62 in our computations. 



3 DISCUSSIONS 

Let us begin by plotting the temporal evolution in 61, 
governed by Eqn. ([l}, that last for couple of years. For 
Fig.[T] we adopted the likely system parameters, a rising from 
the radio and infrared observations (|Champion et al. 2008 ). 
Therefore, we let mo = 1.74 MQ,mi = 1.051 MQ,m2 = 
0.9 Mq,6i = 0.436678,01 ~ 0.211 AU and gi = 141.65°. 
Notice that we have employed the values listed in Table. 1 
of Champi on et al. (2008), arising from the 'DDGR' tim- 
ing model l|Tavlor fc Weisberg 198Sil ). The other quantities 
required to obtain 61 (t) are i,a2 and 62 and as they are not 
constrained from any observations, we treat them mostly as 
free parameters in these plots. For a given value of 02/01, 
we estimate 62 such that the inequalities given by Eqs. ([Sjl 
and (U) are satisfied and i is taken to be in the range so 
that the inner binary can experience Kozai resonance. The 
upper panel of Fig. [1] provides 61 (t) plots for four different 
values of 02/01, while keeping i = 100°, and in the lower 
panel, we have ei{t) plots for three canonical values of i, 
while keeping 02/01 = 600. The striking feature of these 
plots is the decrease in the eccentricity of the inner binary 
with time and the decremental change in 61 during a year 
appears at the fifth decimal place. From Eq. (|lb[) . it should 
be noted that the magnitude of dei/dt, while employing 
the above mentioned system parameters, depends only on 

(^a \/l — 6|^ (1 — 9^), where a = 02/01 and (1 — 6^) is 
allowed to take values between 1 and 0.4 for HT experienc- 
ing the Kozai resonance. However, the restriction that any 
chosen value for 62 should satify the inequalities, ((3)l and 



implies that \/1 — 63 j remains almost a constant 

(~ 10^^) for 02/01 500. This is reflected in the appar- 
ent convergence of 61 (t) curves displayed in Fig. [T] These 
arguments allow us to state that it is impossible for dei/dt, 
as deducted from Eqs. ((TJ, to be substantially smaller than 
~ 10~'^''s~'^ by varying 02/01,62,1 and m2 and we depict 
these conclusions pictorially in Fig. (2] Further, we note that 
that whether 61 decreases or increases depends crucially on 
the initial value for gi and as gi = 141.65° at the present 
epoch, 61 must decrease with time. 

We observe that the above estimates for the change in 
61 over a two year period is substantially higher than the 
reported precision of t he orbital eccentricity for the MSP 
(Champi on et al. 2008h . Moreover, the latest timing results 
for the PSR J1903 + 0327 (Freire, priv. commun.) imply that 
6 ~ 10"i®s-i which is very much smaller than our predic- 
tions based on HT configurations for the MSP. Therefore, 
we rule out the possibility that the MSP binary is part of a 
HT configuration experiencing Kozai oscillations. 

Other measurements can be also used to constrain the 
orbital configuration for the MSP binary. A rough estimate 
for the maximum amplitude of the spin period derivative P 
is obtained by treating the inner binary as a single object 
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of mass mo + mi which for ms a wide binary w ith stellar 
mass companion of mass m2 (jjoshi fc Rasio 19971 ). Further, 
we let the projected semi-major axis of this new binary to 
be 600 times that of the inner binary and choose a value 
for 62 consistent with inequalities Q and Q. The associ- 
ated line of sight acceleration leads to an estimate for the 
maximum amplitude of P of the order of 10~^^s/s and this 
is much higher than the reported value P ~ 10~^''s/s of 
Champion et al. (2008). This too suggests that the MSP 
binary is not part of triple system. With the knowledge 
that the MSP binary is not part of a HT, we note that 
the evolution of the argument of periastron of the binary 
can arise only from the third term in Eq. (|la|l . Therefore, 
the reported rate for the advance of periastron can safely 
be explained with the help of general relativity. However, 
the measured periastron advance rate can originate from 
the classical spin-orbit coupling if the binary companion is 
a main-seq uence star having unusually low ro tational period 
< 0.5 day (|Champion et al. 20081 : IWex 19981 ). 

Champion et.al (2008) had also discussed the possibil- 
ity of further optical observations of the main-sequence star 
to probe the triple nature of the system. Here we briefly ex- 
plore the feasibility of observing the optical band radiation 
of the massive white-dwarf companion in the optical image 
which contains the main-sequence star. Using white dwarf 
models by Lohmer et al. (2004), we find that a white dwarf 
of mass ~ 1 Mq and age ~ 1.8 x 10^ years will have an 
effective temperature T^jf ~ 10* K and an absolute visual 
magnitude M„ ~ 13.0. Taking the extinction coefficient in 
the visual band (A„) to be 4.9 and the distance ~ 6.4 kpc, 
employed by Champion et.al (2008), we get an apparent vi- 
sual magnitude m„ ~ 31.93. Using the white-dwarf colors 
mt — m„ and m^ — m^, as given in Bergeron et al. (1995), 
and ratios of the extinction coefficients (At/Av and Ar/A^j) 
from Cardelli et al. (1989), we get an apparent blue magni- 
tude rrib ~ 33.8 and an apparent red magnitude rrir ~ 31. 
These high values of apparent magnitudes imply that it will 
be quite difficult to observe, in the optical band, the com- 
panion white dwarf in the near future, even though it has 
a different spectral characteristic compared to the mainly 
infra-red optical emission from the main sequence star. 

Let us comment on the possible GC origin for the sys- 
tem, motivated by the fact that more than 15% of the MSP 
binaries in GCs are in ecc entric orbits and s ome of them have 
stellar mass companions (|Champion et al. 2008 ). We adapt 
for PSR J 1903-1-0327 the discussions pursued in Bagchi and 
Ray (2009a, 2009b) to theoretically explain the distribution 
of orbital eccentricity of binary radio pulsars in GCs in the 
light of interaction models. In Fig. |3l we plot theoretical 
separations in the eccentricity-orbital period phase space 
for binaries that could have acquired a given eccentricity 
for its corresponding orbital period through fly-by interac- 
tions in the past 1 Gyr of its existence. Note that these 
separations based on lifetimes depend crucially upon the 
parameter v\o/rn that varies widely between the galactic 
globular cluster systems by over three orders of magnitude 
(here t;io is the velocity dispersion of stars in a GC while 
ni is its central (number) density of stars). The pulsar PSR 
J1903-I-0327 like systems could have been produced in GC 
systems upto a mean i'io/n4 = 0.13 (e.g. a system like NGC 
6752 or even better in Ter 5) but not in any globular cluster 
with a higher value of this parameter. Thus the putative ori- 
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Figure 1. Temporal evolution of the eccentricity of the inner 
binary over two years for system parameters listed in the text. In 
the upper panel, we let a2/ai have values 400, 500, 600 and 800 
and the associated values for 62 are 0.970, 0.977, 0.984 and 0.991, 
respectively. The 62 values are chosen so as to satisfy the two 
inequalities, though for the 02/01 = 800 case, the conservative 
inequality for 02/01, given by Eq. (|4]l, is not satisfied. Note that 
the decrements in 61 are fairly independent of 02/01 for 02/01 > 
500 and this may be attributed to the 1/(1 — e^) dependency of 
the right hand side of Eq. l llbl l. In the lower panel, the mutual 
inclination angle between the two binary orbits i takes values 
45°, 90° and 140° and we fix 02/01 = 600 which leads to 62 being 
0.984. The higher decrement for the i = 90° case is expected due 
to more efficient occurring of the Kozai oscillations, for which i 
should be in the range ~ 39°-141°. The system parameters in 
this panel are the same as those for the upper panel. 

gin of PSR J 1903-1-0327 inside a GC could have taken place 
only in a restricted set of GCs. This further restricts the en- 
vironment for production of such pulsars and therefore the 
overall probability for origin of this pulsar in a GC must 
be smaller than that estimated by Champion et al (2008). 
Note that the fly-by interactions as well as the exchange and 
merger interactions scale with the parameter vio/n4,. 



4 CONCLUSIONS 

We probe the observational implications if the recently dis- 
covered highly eccentric galactic binary MSP PSR J1903 -I- 
0327 is part of a HT experiencing the Kozai oscillations. We 
model the binary MSP and the main sequence star, observed 
in the fleld of the MSP in infrared, to be part of a HT, while 
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Figure 2. Colour coded contours of ei = de\/dt in /sec aris- 
ing from Eqs. JTJ. In the upper panel, we vary m2 and 62 for 
constant a and i values, while in the lower panel a and 62 are 
varied with fixed m2 and i values. The coloured regions are de- 
fined by the inequalities (|3]l and (|4]l, providing the left and the 
right boundaries of the allowed regions, respectively. The chosen 
value of i ensures minimum possible ei, which is much larger than 
the recently measured e value for the MSP binary. 

incorporating the effects due to the quadrupolar interactions 
between the binary orbits and the dominant contributions 
to the general relativistic periastron precession in the inner 
binary (|Blaes et al. 20o3 ). Using resuhs from the timing of 
the MSP, we demonstrate that the eccentricity of the binary 
pulsar should decrease at the present epoch and estimate the 
rate to be ~ 10~^^s~^. As this is clearly inconsistent with 
the measured e ~ 10~'^®s~^ for the binary (Freire, priv. com- 
mun.), we rule out Kozai resonance as the reason for the high 
orbital eccentricity for the binary MSP PSR J1903-)-0327. 
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Figure 3. Fly-by isochrones for tjiy = lO" years in the eccentric- 
ity - orbital period plane for difl'erent values of the GC parameter 
''lo/"4; where i;io is the velocity dispersion of stars in a GC and 
n4 being the number (central) density of stars. The GC binary 
pulsars with spin periods < 20 ms have been plotted with black 
o symbols and the eccentric galactic disk millisecond pulsar PSR 
J1903-I-0327 has been marked with a red o. The region to the lower 
right of each of these lines have shorter interaction timescales. 
This indicates that our binary pulsar could have originated only 
in a rather special GC environment with characteristics belonging 
to the upper left portion of the isochrone with viq/uh = 0.134. 
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